Orbital ferromagnetism and anomalous Hall effect in antiferromagnets on distorted 

fee lattice 
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The Berry phase due to the spin wavefunction gives rise to the orbital ferromagnetism and anoma- 
lous Hall effect in the non-coplanar antiferromagnetic ordered state on face-centered-cubic (fee) 
lattice once the crystal is distorted perpendicular to (1,1)1) or (1,1,0)- plane. The relevance to the 
real systems 7-FeMn and NiS2 is also discussed. 
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It has been recognized for a long term that the chi- 
rality plays important roles in the physics of frustrated 
spin systems These de grees of freedom are dis- 

tinct from the (staggered) magnetization, and could show 
phase transition without magnetic ordering [0-^. Espe- 
cially since the discovery of the high-Tc cuprates, the 
scalar spin chirality 



Xijk — Si ■ {Sj X Sk) 



(1) 



has been a key theoretical concept in the physics of 
strongly correlated electronic systems This spin 

chirality acts as the gauge flux for the charge carriers 
moving in the background of the fluctuating spins. In 
order for the spin chirality Xijk to be ordered, both the 
time-reversal (T) and parity (P) symmetries must be bro- 
ken. Broken T and P symmetries in 2D bring about 
many intriguing physics such as parity anomaly j^J^ , 
anyon superconductivity Q, and quantized Hall effect 
without external magnetic field | [Tc| ]. A physical realiza- 
tion of the last one has been discussed in the context 
of anomalous Hall effect (AHE) in ferromagnets via the 
spin chirality mechanism p^-p^. 

In this paper we explore the chiral spin state in the 
ordered antiferromagnet (AF) on the three-dimensional 
face-centered-cubic (fee) lattice. The AF on the fee lat- 
tice is a typical frustrated system, and nontrivial spin 
structure with the finite spin chirality in eq.(l) is ex- 
pected. For example, in the charge transfer (CT) in- 
sulator NiS2 ^ and in the metallic alloy 7-FeMn 
the non-coplanar spin structure (so-called triple-Q struc- 
ture shown in Fig. la) has been observed. A theoreti- 
cal explanation for this structure is the following. Let 
us consider the case where the lattice points are di- 
vided into 4-sublattices as shown in Fig. la. Denot- 
ing the (classical) spin moment at each sublattice as 
Sa (a = 1,2,3,4), the 2-spin exchange interaction en- 
ergy is written as H2 oc (X]a=i 4 "^a)^- Therefore the 
condition of the lowest energy J2a=i i^a — does not 
determine the spin structure and leaves many degener- 
ate lowest energy configurations. Then the interactions 
which lift this degeneracy such as the 4-spin exchange 



interaction become important |18| , p^ . In paticular the 
phenomenological Ginzburg-Landau theory for the 4-spin 
exchange interaction is given as H4 — J4 J2a^b(^<^ ' ^b)'^ 
pOf . With positive J4, the ground state configuration 

is given by = (tS'TI'TI)' = (tS' "75 ' "Tf)' 

= (-71' 71' -73)' - (-73' -73 '73)' -^-^^ 
each direction corresponds to the four corners from the 

center of a tetrahedron. This non-coplanar spin structure 
gives the scalar chirality in eq.(l) locally. For the itiner- 
ant system 7-FeMn, a recent band structure calculation 
pof concluded the stability of the triple-Q spin strucure, 
in agreement with experiment p7[. 
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FIG. 1. (a) Triple-Q spin structure on fee lattice, (b) Re- 
lation of the 4-spin moments Sa {a = 1,2,3,4). (c) Triangle 
lattice as a cross section of he fee lattice perpendicular to 
[l,l,l]-direction. 

Experimentally anomalous behaviors in the fee AF are 
often observed. For example there occurs mysterious 
weak ferromagnetism (WF) in NiS2 below the second AF 
transition temperature Tn2 The Hall effect in this 

material is also large and strongly temperature depen- 
dent ||2^. In Co(Sa;Sei_j;)2, the AHE is enhanced in the 
intermediate x region, where the nontrivial magnetism is 
realized |^^, and CeSb shows the largest Faraday rota- 
tion ever observed |24[| . 
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As shown below, the triple-Q AF state (Fig. la) on 
the distorted fee lattice gives the orbital ferromagnetism 
accompanied by the AHE and even the realization of the 
three-dimensional quantum Hall liquid ]25| ] without the 
external magnetic field. 

The Hubbard on-site repulsive interaction U can be 
decomposed in terms of the Stratonovich-Hubbard trans- 
formation as 



U -.2 

2^^ 



Uipi 



(2) 



where a are the Pauli matrices and tpi is the spin fluc- 
tuation field. This field acts as the local effective mag- 
netic field for the electron , which can be regarded as the 
c-number in the mean field approximation for the mag- 
netically ordered state. The symmetry properties of the 
ground state is captured by this approximation , and the 
quantized Hall effect discussed below is stable against 
the fluctuation of the (^-field as long as the gap does not 
collapse. Therefore we can consider this system as the 
free electron model in the background of the static field 
Uipi. In paticular when this field U\(pi\ is much stronger 
than the transfer integrals tij^s, we can map the model 
into that of spinless fermions H — ^if fl fj '^i^h 

the effective transfer integral 

= t^^.(cos I cos I -t- sin I sin |e-^(*'-^^)), (3) 

where \xi) = [cos sin ^e*"^*]* is the spin wavefunction 
and 9i and are the polar coordinates of the spin direc- 
tion. The phase factor appearing in f^^ can be regarded 
as the gauge vector potential (r ) , and the correspond- 
ing gauge flux is related to Xijk in eq.(|l|) However, 
it should be noted here that the effective magnetic flux 
vanishes when averaged over the unit cell because of the 
periodicity of the vector potential a^(r) = afj_{r + Ti) (Ti : 
primitive vector). Therefore, the net effect of the gauge 
flux comes from the fact that there are more than two 
atoms and/or orbitals in the unit cell and the resultant 
multiband structure. Then each band is characterized by 
the Chern number [26| . The Chern number appears as a 
result of the spin-orbit interaction and/or the spin chiral- 
ity in ferromagnets 0,^^ . In ferromagnets the T-broken 
symmetry is manifest, while in AF the T-operation com- 
bined with the translation operation often constitutes the 
unbroken symmetry. In the latter case, the nonzero Hall 
conductivity axy is forbidden. However when there are 
more than two sublattices and the spin structure is non- 
collinear, this combined symmetry would be absent and 
finite axy is not forbidden. 

In order to facilitate the understanding, first consider 
the triangular lattice, which is the (l,l,l)-cross section 
of the fee lattice (Fig. Ic), and tight binding model 
with nearest neighbor hopping given in eq.(j3|) has the 



gauge fiux of 7r/2 penetrating each triangle. The unit 
cell is composed of eight triangles including four atoms. 
Thus the total gauge fiux penetrating the unit cell is 
An = 0(Mod(27r)), which is consisitent with the peri- 
odicity of the lattice. The Chern number can be esti- 
mated analytically as — x(7r) lower (higher) en- 
ergy band. Therefore this two-dimensional spin configu- 
ration on the triangular lattice offers an example where 
the spin chirality orders without the ferromagnetic spin 
moment. 

When one considers the three-dimensional fee lattice, 
on the other hand, there are three other cross sec- 
tions equivalent to [l,l,l]-direction, namely [-1,-1,1],[- 
1,1,-1], [1,-1,-1]. Therefore it is naturally expected that 
the net spin chirality is zero, because the spin chiralities 
are the vector quantities and the sum of these four vec- 
tors is zero. Actually a^y = Cyz = <Jzx = for the fee 
lattice as shown below. However, this means that the 
chirality remains finite when the symmetry between the 
four directions is violated. For example, when the lat- 
tice is distorted along [l,l,l]-direction, it is expected that 
(JH ~ cTxy = cTyz ~ Ozx bccomes finite. We express the 
distortion along [l,l,l]-direction by putting the transfer 
integral within the (l,l,l)-plane as iintra — 1, while that 
between the planes as iintcr — \~d [^ . As the unit cell 
is cubic shown in Fig. la , the first Brillouin zone (BZ) 
is cubic : [— -, -]^. From now on, we set a = 1 . Then 
the Hamiltonian matrix i?(A;) for each k is given by 



( e-'6/2 e^e/i h \ 
V h e-*^/i e^^/2 / 



(4) 



where A = 2(1 - d)cos(^ + ^) + 2cos(--^ + ^), 
h = 2(1 - d)cos(^ + ^) + 2cos(-^ + ^), h = 
2(1 - d)cos(^ + \) + 2cos(-^ -I- \). In this Hamil- 
tonian, two bands (£\^2 = —Vfi + /2 + /a) and upper 
two bands (£3,4 = V/i + /2 + /s) are completely de- 
generate. At d — 0, these two dispersions touch along 
the edge of the 1st BZ, i.e. {kx = ±7r, fcy = ±n,kz), 
{kx,ky = ±7r, kz = ±7r), [kx = ±7r, ky, kz = ±7r). Fixing 
kz^ for example, the /c-point {kx — ±TT,ky — ±tt) is the 
center of the massless Dirac fermion (Weyl fermion) in 
(2-f 1)D. Therefore the band touching along the edge can 
be regarded as an assemble of the (2+l)D Weyl fermions. 
What happens for finite d differs for positive and nega- 
tive values of d. For d > (elongation along [1,1,1]- 
direction), all the Weyl fermions along the edges open a 
gap and turn into the massive Dirac fermions. Therefore 
the gap opens in the density of states centered at e = 0. 
For d < (suppression along [l,l,l]-direction), on the 
other hand, all the (2-|-l)D Weyl fermions along the edges 
open the gap; two (3+l)D Weyl fermions emerge instead 

at {kx,ky,kz) = ±2arcsin(./|^)(l,l,l) = ±-0(1,1,1) 
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where ± correspond to right- and left-handed chirahty 
1^ . Here we consider the mass of the Dirac fermions at 
the edges of the BZ. For d > 0, the mass of the (2+l)D 
fermion along the fc^-axis is positive for every fc^, while 
that changes sign at fc^ = ±D for < 0. 

Based on this observation, the transverse conductivity 
can be calculated in terms of the formula in [E6| 



Ef dk^ 
/ ^ 



/(e„(A:))(V^ X i*„(^)). 



2tt 



(5) 



for all while (Jxy{kz) 



O'xyikz) — 

kz = {-D,D) and a^yikz 
and [D.-k] in the case of nagative d < 0. 
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where /(e) is the Fermi distribution function, and 
An{k) = {nk\Wj:\nk) is the vector potential created by 
the Bloch wavefunction \nk). Let us consider fi = 
case, i.e. quarter filled case. axy{kz) is the Hall con- 
ductance in (2-|-l)D system with fixed kz, which is de- 
termined by the parity anomaly caused by the Dirac 
fermions along the edge of the BZ. Explicit calucula- 
tion shows at T = OK , a^yikz) = -^sign{f2{kz) ■ d) 
where /2(fc^)=(4-2d)cos(%)2-h2d-sin(%)2. For d > 0, 



Integrating 



over kz, we obtain cr^y — for d > and (Jxy ~ 



for d < 0, while dxy vanish for 



— arcsinf 

0. One can easily check that Uxy — <Jyz = <^zx — (^Hi 
and that the sign of an changes when all the spins Sa 
are inverted {Sa —Sa)- One can confirm analytically 
in eq. that the cancellation between the contributions 
from the two fc-points related by the parity (P) symme- 
try, e.g. {kx,ky,kz) and, {-kx,ky,kz), occurs and axy 
is zero for the undistorted fee lattice. This remains true 
even when we modify Sa from those given in Fig. lb. For 
example even with finite magnetization {J2t=i 7^ 0), 
axy{uj) = for any uj on the undistorted fee lattice. In 
short, P-symmetry is not broken while T-symmetry is 
broken in the undirtorted fee lattice, where finite an 
is forbidden. An distortion along [l,l,l]-direction does 
break this P-symmetry and produces finite an- 

Next we turn to the half filled case. When the chemical 
potential is in the Mott gap, i.e. Mott insulator , the 
d.c. an vanishes. However the optical aui^o) can be 
finite, which corresponds to the Faraday and/or the Kerr 
rotation. In order to calculate (T//(a;), we have to take 
into account optical transitions between the lower and 
the upper Hubbard bands. In the mean field picture, this 
corresponds to the two splitted bands due to the effective 
magnetic field U(fi in eq.(2). Diagonalizing the tight- 
binding Hamiltonian with this local "magnetic field" , we 
can calculate cr//(a;). (Th{(jj) shown in Fig. ^ is for d — 
~0.1,U\ip,\ = 5. 
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FIG. 2. crH(a;) for d= -0.1,C/|<^i| =5 

Now we turn to the orbital ferromagnetism induced by 
the distortion. Orbital magnetization M is given by 



M = lim 





dB 



dfi. 



(6) 



I.e. 



'xy\ii 



In terms of the Stfeda formula |2£ 
^^yU + ec^^^^l^ with cr^yl^ i^TT[JxG+{ii)Jy5{ii - 
H) — JxS{fJ. — H)JyG~{fj,)], we estimate the integrand in 



eq.(§) from axylfj^ and ct^^I^. 




FIG. 3. Magnetization M as a 
Ulifiil = 5 {we set ;^~iintra = 1 ) 



function of d, where 



Shown in Fig. 3 is the orbital magnetization M as a func- 
tion of the distortion d for the half filled case. The sign of 
orbital magnetization M changes when all the spins are 
inverted. (Magnetization depicted in Fig. 3 corresponds 
to the triple-Q structure in Fig. la,b.) Note that M 
is finite even though the d.c. an = in this case, be- 
cause the former is determined by the integral over the 
occupied states. 

In real materials, the configuration with all the spins 
being inverted (which means that the chirality is also in- 
verted) has the same energy. Thus it is expected that the 
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domain structure of these two chiralities is formed. In or- 
der to align this chiral domain, we must cool down into 
AF order phase by applying external strain and external 
magnetic field, which couples to the orbital magnetiza- 
tion and prefers one of the chiral domain. (For d > 0, 
the magnetic field in [l,l,l]-direction prefers the triple-Q 
structure as in Fig. la,b.) Lastly it is noted here that 
the distortion along [1,1,0] or other three equivalent di- 
rections produces the similar effect as discussed above, 
while that along [1,0,0] etc. does not. This also offers a 
way to test this idea experimentally. 

Now we discuss the experimental relevance of these re- 
sults. One of the most promising materials for this spin 
chirality mechanism is the itinerant AF 7-Fea;Mni_2; al- 
loy as mentioned above [|l^ . In this material, the triple-Q 
structure is observed for 0.35 < x < 0.8. This mate- 
rial remains metallic even below Tjv. Because the crystal 
structure is the undistorted fee, one needs to apply uniax- 
ial stress toward the [1,1,1] or [l,l,0]-direction. Although 
the band structure is rather complicated and there is no 
gap in the density of states [^ , we expect the transverse 
conductivity an of the order of e^d/ha = lA00d^l~^cm~^ 
for a = 3.6A when the chiral domains are aligned by the 
field cooling. 

Another candidate is the CT insulator AF NiS2. In 
this material the valence of Ni is 2-1- , and the angular 
orbital moment is quenched with d^ configuration. Thus 
spin-orbit interaction is expected to be small, and this 
material is an ideal laboratory to study the spin chirality 
mechanism. It shows two successive magnetic phase tran- 
sitions at T/vi = 40K and Tjv2 — 30K. Between Tjvi and 
T/V2 7 the magnetic strucure is given by Fig. la,b (type 
I AF state). At the type II AF structure appears 
in addition to the type I structure. This is accompanied 
by the rhombohedral distortion and with the mysterious 
WF [H. This lattice distortion {d > 0) lifts the quasi- 
degeneracy of type I and type II structures jso). This 
brings about finite orbital ferromagnctism, which can be 
detected by hysteresis in the magnetization curve under 
the magnetic field. Therefore the present results give a 
possible scenario for the WF in NiS2. Even in the tem- 
perature range between Tjvi and Tpf2, where lattice is 
not distorted spontaneously , suppression toward [1,1,1] 
or [l,l,0]-direction is predicted to bring about an orbital 
ferromagnctism. (The compressibility of this material is 
of the order of lO'^kbar"^ and we need ~140kbar 
uniaxial pressure to produce d = —0.1 for example.) Al- 
though the d.c. ani^j) is zero at T = OK, the optical 
cth(w) is expected to be finite, as shown in Fig. 2. As for 
other fee AF such as Co(Sa;Sei_a;)2, detailed studies on 
the spin structure by neutron scattering are highly desir- 
able, with which the mechanism of the Hall effect can be 
dictated. 
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